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Abstract 

Ten years of operations of the Swift satellite have allow us to collect a small sample of long Gamma-Ray Bursts 
(GRBs) at redshift larger than six. I will review here the present status of this research field and discuss the possible 
use of GRBs as a fundamental new tool to explore the early Universe, complementary to quasar and galaxy surveys. 
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1. Introduction 

In the standard ACDM scenario, the first stars, the 
so-called Population III (PopIII) stars, are predicted to 
form in dark matter minihalos of typical mass ~ 10® Mq 
at z ~ 20 - 30 out of a gas of pristine composition 
(see Ciardi & Ferrara 2005 : Bromm & Yoshida 2011 


for recent reviews). Thanks to their peculiar chemical 
composition, they are expected to be more massive than 
the subseq uent stellar populations with typical mass of 
~ 40 Mq ( Hosokawa et al.L 201 ih. possibly exte nding 
up to hundred solar masses dHirano et al. . 20141) . The 
formation of these stars mark the fundamental transi¬ 
tion from a simple, very homogeneous Universe to the 
complex and structured one we see already in place a 
billion years after the Big Bang. During this period of 
time two fundamental transitions are expected to occur: 
(i) a change in the SFR mode, i.e. from massive PopIII 
to solar-size PopII/I stars and (ii) the cosmic reioniza¬ 
tion, i.e. the change of the inter-galactic medium (IGM) 
from a neutral to a fully ionized state. There is a general 
consensus that the first tra nsition is driven by th e so- 
called chemical feedback ( Schneider et all 2002h . i.e. 
the enrichment of star forming clouds by the first su¬ 
pernova e xplosions above a criti cal threshold of Zcut = 
10®** ( Schneider et aO. 2003 ). As the chemical feed¬ 

back is essentially a local effect, on cosmological scale 
the two populations should be coeval for a long period 
of time (Schneider et al.L 2006 : Tornatore et al. , HoO? 


Maio et akT 2010 ). Understanding how and when the 


change in the SFR mode happened is one of the main 
goals for current studies of galaxy formation in the early 


Universe and the detection of PopIII stars one of the 
main challenges of the next generation of sp ace and 


mam cnallenges ot me next generation ot sp ace anc 
ground facilities (see Bromm & Yoshidal 201 ll for a re 
view). The cosmic reionization of the IGM has been ex¬ 
tensively studied in the last years bo th from a theoretical 


and o bservational point of view (see lLoeb & Furlanetto 


20 nl for a recent review). While many steps towards 
our understanding of reionization process have been 
done, still many fundamental details are missed: In 
which way does it proceed? How gradual and how pro¬ 
longed was the process? Was radiation from early stars 
sufficient to sustain this phase transition? Do PopIII 
stars or quasars have a major role in driving the process? 
Is some other more exotic process at work? etc. 

Historically, the exploration of the distant Universe 
has been carried out following two main pathways: the 
observations of bright q uasars detected in wide shal¬ 
low surveys ( Fanl 2012h . and of distant galaxies iden- 
tified through the dro p-out technique in small fields 
( Bouwens et al.L 20141) . Thanks to their extreme bright¬ 
ness Gamma-Ray Bursts (GRBs) represent an alterna¬ 
tive way to access those early e pochs. As demonstrated 


by GRB 090423 at z = 8.2 (ISalvaterra et all 120091 


Tanvir et ^ 20091) . they can be detected even at dis¬ 
tances much larger than any other cosmic object. In 
principle, their afterglow emission can be observ ed up 
to z ~ 20 ( Ciardi & Loebl 2000l: Gou et al.[ 2004) pro¬ 
viding useful information about the ionization and metal 
enrichment history of the early Universe. Here, I will 
review the present observational status of this research 
field and discuss the possible role of GRBs in the explo- 
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ration of the Universe during and before the reionization 
epoch. 

This paper is organized as follows. In Section 2, a 
brief summary of the observations of the GRBs at z > 6 
and of their host galaxies is given. Section 3 presents 
the expected rate of high-z GRBs both from PopII and 
PopIII stars and the expected properties of the relative 
host galaxies. In Section 4, the use of GRBs as a probe 
of the early Universe is reviewed. Finally in Section 5 I 
present some ideas for the future of this research field. 
The conclusions are drawn in Section 6 . 


2. Observations 


2.7. The Swift high-z GRB sample 

In ten years of operations Swift has detected a hand¬ 
ful of bursts with spectroscopically confirmed redshift 
larger than 6 . In addition, other three GRBs have well 
constrained photometric redshift above this limit. The 
observed high-z sample represents ~ 1% of all Swift 
bursts, ~ 2.5% of those with known z. The main prop¬ 
erties of the z > 6 GRB sample are given in Table 1. 


GRB 050904 at z=6.3 Kawai et al. ( 20061) 


This burst was firstly imagined by the 25-cm tele¬ 


scope TAROT (IKlotz et al.L 120051) . Its high-z na 


ture was rec ognized by multi-wavelength photo 


metri c data ( Haislip et al. . 20061: Tagliaferri et al 


2005h and firmly confirmed spectrocopically three 
days a fter the Swif t trigg er by the Subaru tele¬ 


scope (IKawai et al.L l2006h . The afterglow spec¬ 


trum provided a upper limit on the neutral hydro¬ 
gen fraction at th e GRB redshift of xhi < 0.17 


at Icr confidence (iTotani et al.L 120061) and a mea- 
sure of the metalli city at the level of ~ 0.1 


Z„ 


(IKawai et al.Ll2006fl . Recently. iThone et al.l (120131 


revised this value, inferring a slightly lower metal¬ 
licity from the S II 41243 equivalent width, result¬ 
ing in log(Z/Z 0 ) = -1.6 + 0.3. The afterglow spec¬ 
tral energy distribution (SED) requires the pres¬ 
ence of SMC or supernov a (SN) type dust a t a 
ley el of Ay = 0.15 + 0.07 ( Stratta et al. 2011 but 
Zafar et al1l201 Ibh . The modeling of afterglow 


see 


data from X-ray to radio suggests GRB 050904 
to be an energetic burst blowi ng up in a dense 


medium with n ^ 680 cm ^ (iFrail et al.l 12006 


Laskar et al.L 120141). 


GRB 080913 at z = 6.7 ( Greiner et al. . 20091) 


The high-z nature of this burst was recognized yia 
the detection of a spectral break between the i' 


and z' bands of the GROND instrument and then 
confirmed spectroscopically by VLT obseryations. 
The analysis of the red damping wing constrained 
+ht < 0.73 at 90% confidence leyel ( Patel et al 


2010l) . In spite of the rapid follow-up campaign, 
the faintness of the afterglow preyented the detec¬ 
tion of any m etal absorption lin e, but the Su+Siu 
at 2.9cr leyel ( Patel et ^ 2010t) . Dust absorption 
of A y - 0.12 ± 0 .03 is found from SED fitting 


(IZafar et al.l 1201 lah . 


GRB 090423 at z = 8.2 (ISalyaterra et all 12009 
Tanyir et al. . 2009h 


The spectrosc opic redshift of this b urst was se- 
cured by TNG ( Salyaterra et al. , 2009 ) and by VLT 
( Tanyir et al. . 20091) . and still represents the dis¬ 
tance record for a cosmic object. Radio obser¬ 


vations by VLA were reported by IChandra et al.l 
( 2010h . The analysis of the multi-wavelength 
dataset show that the afterglow is reminiscent of 
many other lower redshift bursts, suggesting that 
in spite of its extreme redshift, its progenitor and 
the medium in which it Flowed up were not pe¬ 
culiar. Indeed, its detection is consistent with 
the h igh-z tail of PopII/I G RB redshift distribu¬ 
tion ( Salvaterra et al.L 20091) . From the X-ray to 
optical S ED no absorption by dust is evident with 
Ay < 0.1 (Salvaterra et al.l20()9 : Tanvir et al.l2009 


but see Laskar et al.ll2014l) 


Castro-Tirado et al 


= 5.9 (Chornock et al.. 

2013 

2013; iTotani et al.L 

2014 


Hartoog et al.Ll2014h 

The redshift o f this 
Gemini-North (Chornock et al.l 120131). 


burst 


was 


obta ined by 
GTC 


( Castro-Tirado et^kT 2013 b. Subaru dTotani et al 


20141) and VLT ( Hartoog et al. . 2014 ). In par¬ 


ticular, the superior resolution and wavelength 
coverage of the VLT/X-shooter instrument showed 
the potentiality of GRBs as tool to study in great 
details the metal enrichmen t of star formin g regio n 


inside high-z galaxies (IHartoog et al.L gOlAl). 


Precise column densities of H, Al, Si and Fe are 
reported together with limit on C, O, S and Ni. 
The host metallicity is constrained to be in the 
range of 0.03-0.06 solar and the high [Si/Fe] in 
the host suggests the presence of dust depletion 
(though Ay < 0.05 from SED fitting). The best 
fit of the Lya absorption line is obtained for 
log(VHi) = 19.94 + 0.01 and negligible neutral 
hydrogen in the external medium, with xhi < 0.03 
at 3cr significance. 
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GRB 

z 

Ep 

F 

^ISO 

log(AHl) 

log(AH,x) 

Z 

Ay 

^UV,host 

SERhost 



[erg] 

[erg] 

[cm“^] 

[10^' cm“^] 

[Zo] 


[AB] 

[M©yr->] 

050904 

6.3 

3178 

1.24 X 10^"^ 

21.6 

63!g 

-1.6+ 0.3 

0.15 + 0.07 

> -19.95 

<4.1 

080913 

6.7 

1008 

7 X 10^2 

19.84 

95+89 

10^ 

- 

0.12 + 0.03 

> -19.00 

< 1.3 

090423 

8.2 

746 

1.88 X 10^3 

- 

- 

<0.1 

> -16.95 

<0.38 

130606A 

5.9 

2028 

2.7 X 10^^ 

19.93 

< 30 

-1.35 + 0.15 

<0.05 

- 

- 

140515A 

6.3 

376 

5.1 X 10^2 

18.62 

< 226 

< -0.8 

0.11+0.02 

- 

- 

090429B 

9.4 

437 

4.31 X 10^^ 

- 

140+10 

- 

0.10 + 0.02 

> -19.65 

<2.4 

120521C 

6.0 

- 

1.9 X 10^3 

- 

<60 

- 

<0.05 

- 

- 

120923A 

8.5 

376 

5.1 X 10^2 

- 

< 720 

- 

- 

- 

- 


Table 1: List of the GRBs at z > 6 detected by Swift and of their observed properties: peak energy (Ep), isotropic equivalent energy in the 1-10000 
keV range (£'iso). hydrogen column density in the host (A^hi)^ X-ray equivalent hydrogen column density (A^h.x )5 host metallicity (Z) and dust 
extinction (Ay). The last two columns report the limits on the host galaxy luminosity (Muv,host) ^nd SFR. See text for references. 


GRB 140515A (iChornock et al. 


2014 


Melandri et al.Ll2015h 


The redshift of this burst has 


by Gemi ni-North dChornock et al. 


and VLT (IMelandri et al. 


201 5h . 


been secured 

2014. GTC 


Chornock et alJ 


(120141) analysed the Gemini-North spectra finding 
no evidence of narrow absorption lines, indicat- 
ing a host me t allicity Z < 0.15 Z©. However, 


Melandri et al.l ( 12015h by modeling the X-ray to 


optical SED found evidence for dust absorption to 
the level of Ay ~ 0.1 indicating some metal enrich¬ 
ment. The red damping wing of Lyman-a can be 
fitted equally well by a single host galaxy absorber 
with log(NHi) = 18.62 + 0.08 or a pure IGM ab- 
sorption with neutral h ydrogen fraction jchi ~ 0.06 
( Chornock et al.L 2014 ). 


well fitted assuming constant-density medium with 
no ~ 0.05 cm“^ and no dust extinction. Ay < 0.05. 
The radio observations revealed the existence of a 
jet break corresponding to an jet angle fljet - 3 de¬ 
grees. An isotropic energy Eiso = 1.9 + 0.8 x 10^^ 
erg is obtained by extrapolating the BAT fluence. 


GRB 120923A at z ~ 8.5 (iTanviiibOlsI) 


The optic al-NIR afterglow of this burst was imag¬ 
ined by Gemini-North telescope ~ 1.4 hours af¬ 
ter the Swift trigger. The blue IR color, H-K~ 0.1 
mag together with the break between the Y- and 
J-band is suggestiv e of a strong high-z candidate 
( Levan et dl 20121). The fit to all photometric data 


gives z ~ 8.5 (iTanviri 120131) . No radio or mm de¬ 
tection is reported. 


Other three GRBs have accurate photometric redshift 
measurement that place in the z > 6 sample: 


GRB 090429B at z ~ 9.4 (ICucchiara et al.L 120111) 


Cucchiara et al.l (1201 ll) collected the afterglow data 
obtained with Gemini-North, VLT and GROND. In 
the best fit model these data are all consistent with 
a photometric redshift of z = 9.4 and low extinc¬ 
tion Ay = 0.10+0.02. A secondary solution at very 
low redshift is still allowed by the SED fitting, but 
it seems unlikely due to the lack of any detection 
of the GRB host galaxy (see next Section). 


GRB 120521C at z ~ 6 (iLaskar et al.L 12014 


The large multi-wavelength dataset of the after¬ 
glows lead to z - 6.0 also supported by the analysis 
of the very low signal-to-noise Gemini-North spec¬ 
trum. All afterglow data after ~ 0.25 days can be 


2.2. General properties ofhigh-z GRBs 

In spite of the burst to burst differences and of the 
small size of the sample, some interesting clues about 
the nature of high-z bursts can be gathered by comparing 
them with bursts at lower redshifts. Thus, in the follow¬ 
ing, I will compare the properties of z > 6 bursts with 
those of a well selected, complete sam ple of bright Swift 
bursts (BAT6, Salvaterra et al. . 2012h . The BAT6 sam¬ 
ple has been selected only on the basis of the y-ray peak 
flux as seen by 5w//f/BAT and is charac terized by a very 


high completeness in red shift (95%, ISalvaterra et al 


2012tlCovino et al.L 120131) . Thus, it represents a perfect 
sub-sample to compare the results of the high-z burst 
population with. 

The top-left panel of Eig. 1 shows the position of 
z > 6 bursts in the Ep - E^ c, plane. The blac k points 


represents the BAT6 sample (iNava et al.Ll2012h and the 
shaded area shows the 3cr scatter around the - Ek„ 
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correlation. High-z bursts nicely follow the Amati cor¬ 
relation suggesting that their prompt emission proper¬ 
ties do not differ significantly from bright lower-z bursts. 
This consistency calls for a common central engine (and 
progenitor star). 

In the top-right panel of Fig. 1 1 have then compared 
the rest-frame 2-10 keV light curves normalized for the 
corresponding Eiso of z > 6 bursts with those obtained 


from the BAT6 sample (iD’Avanzo et al.Ll2012h . Indeed, 


the X-ray light curves of low-z GRBs tend to cluster 
when such a normalization is performed and strong cor¬ 
relation between Lx and E^so is found at early times. 
This effect is seen also in high-z GRBs and their X- 
ray light curves can not be distinguish from those of 
low- and intermediate-z bursts. This conclusion is also 
supported by more extensive, although model depen- 
dent, analysis of afterglow emission. It ha s been shown 
( Chandra et al.L 2010l: Laskar et all 2014 ) that the total 
energy, microphysical parameters and medium proper¬ 
ties of the high-z GRBs are in line with those derived 
for lower redshift bursts. 

In the bottom-left panel of Fig. 1, the available 
host dust extinction estimates for z > 6 bursts are 
plotted against thos e obtained for the BAT6 sample 
( Covino et al. . 2013h . Being selected only on the basis 
of the GRB y-ray peak flux, the BAT6 sample is repre¬ 
sentative of the extinction distribution in GRB selected 
galaxies: 50% of bursts suffer less than 0.3-0.4 mag 
extinction and only 13% of GRBs have Ay > 2 mag. 
The high-z sample is characterized by little or no ab- 
sorptiorQ, suggesting a decrease in dust co ntent in star ¬ 
formin g environments at high redshifts (IZafar et al 


I 2 OI lah . Similarly, the metallicities inferred for z ~ 6 
bursts are in line, though at the lower end, with the dis- 
tributio n of Z measured in l ower redshif t GRB after¬ 
glows ( Savaglio et al.L 2009; Snarre et al. . 2014 ). The 
observed Z and A y values (or limits) are indeed in agree¬ 
ment with those expected for high-z galaxies populat- 
ing the faint end of t he luminosity function (see e.g. 


Salvaterra et al.L 1201 Ih . 


At variance with lower redshift bursts, all high-z 
GRBs have measured (or limit consistent with) very 
high X-ray equivalent hydrogen column densities, Ah,x 
( bottom-right panel of Fig. 1). At lower redshift, a 
high Ah.x is usually connected t o highly extin c ted af 


terglows ( Camnana et al.[ 2012 : Covino et al. . 2013 ). 


However, at higher redshifts this interpretation seems 
to be unrealistic given our knowledge of galaxy evo- 


* However, tlie laclt of liighly extincted bursts may reflect a bias in 
the sample as even a small amount of absorption at z > 6 will prevent 
the detection of the optical afterglow and thus of the redshift measure. 


lution at early times. Moreover, little or no dust ab¬ 
sorption is found in high-z bursts, calling for unusual 
metal-to-dust ratios. In recent years, a well defined 
trend of increasing Ah,x with r edshift has been found 
using both GRBs and AGNs ( Camnanaet^ ^ 20jT ; 


Behar et al. . 2011; Starling et al. . 20131: Camnana et al 


2015h . High-z GRBs extend this relation to higher red¬ 
shifts quite nicely with no A^h.x value or limit belo w the 
observed lower envelope. Camnana et alJ ( 2015 ) have 
shown, by means of dedicated numerical simulations, 
that this effect can be naturally explained by the absorp¬ 
tion of intervening metals along the line-of-sight (LOS), 
being the observed lower envelope due to metals present 
in the IGM. An additional contribution in most of the 
LOSs comes from metals in the diffuse gas in galaxy 
groups at z ~ 0 - 2. 

In conclusion, high-z GRBs do not show any pe¬ 
culiar feature with respect to low- and intermediate-z 
events. Both the energetics and the afterglow properties 
are found to be very similar. The medium in which they 
blow up does not differ too much in terms of density, 
metallicity and dust content. Finally, also the duration 
of their prompt emission is consistent with those of the 
low-z sample takin g into account co s molog ical and in¬ 


strumental effects (iLittleiohns et al.L l2013h . All these 


findings support the idea that the GRBs detected so far 
at z > 6 represent the high redshift tail of the PopII/I 
GRB redshift distribution. 

2.3. Eligh-z GRB host galaxies 

In the last few years, deep s earches with both 
ground- an d space - based faciliti es ( Tanviretall 2012: 


Basa et all 1201 2t [Walter et all 120121 : iBerger et al 


20141) have been carried out searching for the host 
galaxy of high-z GRBs. Indeed, the a-priori, precise 
knowledge of the position and distance of the galaxy 
provided by the GRB detection allows to optimize the 
te kscope setup , pushi ng the instruments to their limit. 

Tanvir et al.l (120 12h reported the observation of the 
field of four GRBs at z > 6 with HST with filters 
choosen to cover the rest-frame wavelength range 1200- 
1500A, i.e. above the Lyman-a. In spite of the very 
deep limits reached (see Table 1), none of the targets 
have been i dentif ied. Similar results were obtained by 
Basa et alJ (120 12h using slightly shallower V LT obser¬ 
vations. More recently, Berger et al. ( 20141) imagined 
the field of GRB 090423 with Spitzer at 3.6 pm and 
ALMA at 222 GHz, but, also in this case, no detection 
is reported. These results imply that all the hosts lie be¬ 
low the characteristic luminosity value at their respec¬ 
tive redshifts, with star formation rates SFR< 4 Mq yr“^ 
in all cases. In particular, GRB 090423 was possible to 
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Figure 1: Comparison of the properties of z > 6 bursts (red squares and lines) with those of a well selected complete sub-sample of bright Swift 
GRB, the BAT6 sample (black points and lines: ISalvaterra et al.l2Q12l) . Top-left panel: correlation where BAT6 bursts are from lNava et al.l 

<2012h . The shaded regions represents the 3(T scatter around the best-fitting relation. Top-right panel: rest-frame 2-10 keV light curves normalized 
to their isotropic energy where BAT6 data are from lD’Avanzo et all (lOl^)- The shaded region represents the 2cr scatter around the mean value 
of the Lx/-£'iso disti 'ibutions at the given rest-frame time. Bottom-left panel: dust extinction as a function of redshift, where BAT6 data are from 
ICovino et^ Bo ttom-r ight panel: intrinsic X-ray equivalent hydrogen column densities Ah,x as a function of redsh ift, where BAT6 data 

are from lCamnana et al.l J2012I) . The shaded regions represent the effect of intervening material along the line-of-sight fsee lCampana et aT]l2015l 
for the details). 
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deriv e a strong limit on t he unobscured SFR< 0.38 Mq 
yr“^ ( Tanvir et al, . 2012 ). while the ALMA non detec- 
tion required the o bscured SFR to be < 5 Mq yr * 
( Bergere^^ , [2^^. 


Tanvir et al.l (120121) also stacked all images deriving 


a limit on the mean SFR per galaxy of < 0.17 Mq 
yr^', consistent with the idea that the bulk of star for- 
mation activity is miss ed in current deep HST surveys 
(ISalvaterra et al.L 1201 ll) . Moreover, these findings offer 
independent evidence that the galaxy luminosity func¬ 
tion is evolving rapidly with redshift with a steeper 
faint-end slope or with a decreasing charac teristic lu¬ 
minosity at high redshift (iTanvir et al.Ll2012h . 


3. Theory 



3.1. The high-z GRB population 

As described in previous Section, Swift has detected 
almost one GRB at z > 6 per year. Considering the 
1.4 sr field of view of the BAT instrument, this cor¬ 
responds to a rat^l of ~ 0.6 bursts yr ' sr '. It is 
worth to ask whether such a rate is expected on the 
basis of our knowledge of the intrin sic GRB redshift 
distribution and lum in osity f unction ( Bromm_&^^oeb 


2007 


Salvaterra et a 

.1 2009bl; Butler et al., 

2010 

Wanderman & Piran 

2 OI 0 I 

iRobertsor 

& Ellis. 

2012 

Salvaterra et al.. 2012 

Ghirlanda et al., 

201. 5b. 


Figure 2 shows the distribution of peak fluxes of 
bursts atz > 6 as measured in the 15-150keV Swift/BAir 
band. The histogram is compared with the model results 
under different assu mptions on the evolution of the GRB 
luminosity function ( Salvaterra et al.L 2012h . The mod¬ 
els have been calibrated by jointly fitting the differen¬ 
tial peak flux distribution of BATSE long GRBs and the 
redshift distribution of the BAT6 sample. The models 
reproduce the rate and peak flux distribution of z > 6 
burst remarkably welj^ with a small excess in the pre¬ 
dicted population at the faintest fluxes. However, since, 
in general, to faint y-ray fluxes correspo nd faint opti¬ 
cal afterglows (e.g. Melandri et akl 2014 ). some high-z 
faint bursts can be missed being the detection of its af¬ 
terglow difficult even with 8-m class telescopes. 


^Wliile tliis represents a lower limit, given the fact that some high-z 
burst can be among GRBs without redshift measurement, it should be 
noted that there is a possible positive bias at work here. Indeed, most 
of the observational programs at large ground- and space-facilities 
have as their primary goal the detection of high-z events. 

^Since the BAT6 sample is limited to z < 5.4 there is no assump¬ 
tion in the model about the peak flux distribution of z > 6 GRBs 


Figure 2: Peak flux distribution of bursts at z > 6 as measured 
in the 15-150 keV Swift/BAT band. The solid (dotted) histogram 
shows the distribution of bursts with spectroscopic (spectroscopic and 

S ihotometric) redshift. The shaded regions report model predictions 
Salvaterra et all l2012h . taking into account the uncertainties in the 
determination of the GRB luminosity function and evolution. Dai'k 
(light) color refers to the density (luminosity) evolution models. The 
dashed line is obtained assuming no evolution and clearly underesti¬ 
mates the number of high-z detections at all peak fluxes. 


By extrapolating model results to faint er fluxes 
( Salvaterra et al. , 2008 ; Ghirlanda et al. , 201.5 ), it can be 
shown that bursts z > 6 represent ~ 10% of the whole 
population, suggesting that GRBs are quite efficient in 
selecting high-z objects. The detection of these events 
is one of the main goals of any future GRB mission. 


3.2. GRBs from PopIII stars 


Different authors 

(W 

eszaros & Rees 

, 201 c 

Komissarov & Barkov 


2010l Suwa & loka 

2011; Tomaetal.. 

201 ll 

Nagakura et a 

ll 2012 

Nakauchi et al., 2012; 

Piro et al., 2014b have nronosed 


that the conditions for jet breakout could be met also 
during the final phases of the collapse of a PopIII 
massive star. Although under different frameworks, all 
models consistently predict PopIII GRBs to be very en¬ 
ergetic events, with total energies exceeding by orders 
of ma gnitude thos e expe cted in PopII events. In par¬ 
ticular Toma et ^ ( 201 ih suggested that PopIII bursts 
could have an equivalent isotropic energy ~ 10^®“^^ 
erg making their detection possible even at the highest 
redshifts. A much longer prompt emission, with typical 


models. However, all these characteristics are shared 
observationally with the population of ultra-long GRBs 
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recen tly detected at much lower redshift (ILevan et al. 
20141) and likely associated with PopII blue supergiant 
stars ( Nakauchi et all 2013h . Thus, they do not repre¬ 
sent a unique feature to firmly indentify a PopIII event. 
Strong evidence could be provided by the absence 
of metal absorption lines in the afterglow spectrum 
down to a level of the critical metallicity. However, the 
measure of such low metallicity (or a limit consistent 
with it) requires extremely high signal-to-noise spectra 
that can be even beyond the capabilities of 30-m class 
telescopes. In the absence of any other observational 
tool to identify PopIII events, radio observations can 
provide an import ant clue. Indeed extremely powerful 


radio afterglows (Ciardi & Loebl 20001 : Toma et al._ 
201 ll: iGhirlanda et all ^13bl) are exp ected giyen the 


peculiar energetics of PopIII events. IGhirlanda et al. 


( 2013bl) showed that PopIII radio afterglows should 


reside in a very distinct position in observed peak 
time-peak flux plane, reaching much larger peak fluxes 
at later times than any PopII GRBs. 

We have seen in Section 12.21 that it is un¬ 
likely that the observed z > 6 bursts are asso¬ 

ciated with PopIII progenitors. This piece of in¬ 


pected rate of PopIII events (Bromm & Loebl 2006 

Came 

si et al. . 201 ll 

Toma et al.. 2011; de Souza et al. 

291 1 ; 

Maio & Barkov. 12014; 

Mesler et al. 

2014). Bv 


means of dedicated numerical sir nulations includ i ng the 
effect of the chemical feedback, Campisi et al. ( 201 ih 


computed the expected redshift distribution of PopIII 
GRBs. The PopIII GRB rate is found to increase from 
z = 20toz = 8-10 and then it decreases rapidly be¬ 
coming negligible at z ~ 5. Under the assumption that 
all PopIII GRBs can be detected by the maxi¬ 

mum allowed rate is ~ 0.03 Gpc“^ yr"'. Assuming a 
typical jet angle of ~ 7 degree and a typical PopIII mass 
of ~ 100 Mq, this corresponds to < 1 PopIII GRB ev¬ 
ery 500 P opIII stars, comparable w ith the PopII/SN type 
Ib/c ratio ( Ghirlandaet al. . 2013al) . 


In conclusion, even if GRBs can really result from 
the collapse of a massive PopIII star, they will be ex¬ 
tremely rare with < 0.05 event yr ' sr * at z > 6. PopIII 
GRBs becom e the dominant observ ed population only 
atz > 10-15 ( Campisi et al. . 2011 ). On the other hand, 
even a single detection of a GRB powered by PopIII 
stars, i.e. a extremely high-z, very energetic, long event 
with a bright radio afterglow peaking a late times, would 


^Tliis may b e not tlie case since the Swift trigger algorithm 
iLien et alll2014f) is not very sensitive to extremely long GRBs. 





Figure 3: Properties of simulated GRB hosts in the redshift range 
Z = 6 — 10. Panels from top to bottom show the distribution of UV 
absolute magnitude, of the SFR (in Mq yr“*), and of gas phase metal- 
licities (in solar units). The dotted line shows the normalized cumula¬ 
tive distributions (see right y-axis). The an'ows i n the first two panel s 
report the limits on Muv and SFR of GRB hosts iTanvir et alll2012ll . 
The shaded areas in the bott om panel correspon d to the metallicity 
measured in the GRB 050904 iThOne et al.ll2()1.3ll and GRB 130606A 
iFlartoog et alil20T^ afterg low spectrum. The arr ow shows the limit 
infen'ed for GRB 140515A IChornock et al.ll2014h . 


represent a breakthrough in the study of the first gener¬ 
ation of stars. 


3.3. High-z GRB host galaxies 

In absence of any detection (see Section 12.31 1. the 
relation of GRB selected galaxies with the typical 
high-z galaxy population has been studi ed by means 


of d edicated semi-analytical models (iTrenti et al 


2 OI 2 I) or numerica l simulations ( Salvaterra et al. , 2013 


Elliott et al.. 2015). 


1 


alvaterra et al 


(120131) derived the expected proper¬ 


ties of the GRB hosts at z = 6 - 10 by means of cos¬ 
mological numerical simulations including all relevant 
feedback effects at these redshifts. GRBs are found 
to explode in bursty galaxies with typical star forma¬ 
tion rates SFR^ 0.03 - 0.3 Mq yr“', stellar masses 
M* ^ 10® - 10^ Mq, specific star formation rates 
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sSFR^ 3-10 Gyr"'. The distribution of their UV lu¬ 
minosities (see top panel of Fig. 3) places them in the 
faint end of the galaxy luminosity function, below the 
capabilities of current instruments. This is consistent 
with the lack of detection and suggests that deep obser¬ 
vations with the James Webb Space Telescope (JWST) 
are required to pinpoint them. Finally, the bulk of 
the GRBs are found to explode in galaxies already en¬ 
riched with metal to Z ^ 0.03 - 0.1 Zq (bottom panel 
of Fig. 3), consistently with the metallicity inferre d in 
z ~ 6 GRB afterglow sp ectra ( S alvaterra et al.ll20 13 . but 


see 


Cen & Kimmll201^ . 


The same simulations allow to predi ct the properties 


of ho st galaxies of PopIII GRB events (ICampisi et al 


201 Ih . These bursts are found to reside typically in ob¬ 


jects at the lower end of the stellar mass with M* < 
10 ^ Mq as they have the higher probability of having 
still a pristine composition. While a metallicity below 
Zcrit is obviously required in the star forming cloud from 
which the PopIII GRB arises, some of the PopIII GRB 
hosts are found in galaxies with Z ~ 10“^ Zq. This 
may reflect the fact that some PopIII event can be asso¬ 
ciated with pock ets of metal-free gas in the outskirt of a 
enriched galaxy (ITornatore et al.ll2007h . 


4. High-z GRBs as a tool 

4.1. Independent measure of the SFR 

It is well established that lon g GRBs are associated to 
the death of massive stars (see IHiorth & Blooml 12012 . 
and references therein). Therefore, they should trace to 
some extent the star formation activity through cosmic 
times and they can be use to measu r e in an independent 
way the global SFR (Kistler et al. . 20091: Ishida et al 


2011 : Robertson & EllisT 2012 ). In this respect, GRBs 


have many advantages compared with usual probes: (i) 
they are detected at higher redshifts; (ii) they are inde¬ 
pendent on the galaxy brightness; (iii) they do not suffer 
of usual biases affecting optical/NIR surveys. Indeed, 
although hampered by small statistics, the detection of 
the few high-z GRBs already suggests that the global 
SFR at z - 8 - 9 is 3-5 times higher than deduced from 
high-z galaxy searc hes through the drop-out technique 


dishida et al.U201 11) . An important caveat to be mention 


here is that GRBs may be not perfect tracers of the SFR. 
Indeed, in order to reproduce t he redshift distribution of 


2007 


2012 

, 


Wanderman & Piran, 2010t Robertson & Ellis 


Salvaterra et al.L 120121) . some kind of evolution 


with redshift is required. Thus, a precise knowledge of 
the nature and the value of this bias is needed to properly 


use GRBs as SFR probe in the early Universe. This can 
be achieved by measuring the GRB LF and its evolu- 


2012: Wanderman & Piran. 2010t Pescalli et al.. 201.51) 

and by studying well selected sam 

pie of low-z GRB 

host galaxies ( 

Boissier et al., 

2013 

Perley et al.l 2013 

Huntetal., 20141 Trend etal.1 12014; Vergani etal. 


20141) . In particular, IVergani et al.l d2014h have shown 
that the stellar mass distribution of a small, but com¬ 
plete, sample of GRB selected galaxies at z < 1 is con¬ 
sistent with the existence of a mild metallicity bias in the 
GRB hosts, with bursts forming more efficiently in en - 
vironments with Z < 0.3 -0.5 Zq ( Vergani et al. . 2014). 
If these findings are confirmed by larger samples, GRBs 
can be considered good tracers of the SFR at least at 
z > 3 - 4. 

4.2. The reionization history 

Similar to quasars, GRBs have been used to mea¬ 
sure the neutral hydrogen fraction in the IGM by fit¬ 
ting the red damping wing of the Lya and therefore 
to place strong constraints on the reionization redshift 


( Totaniet^ , 200^ Chornock et al. , 201^ Totanie]^^, 
2014 ; Hartoog et al. , 2014 ; Chornock et al. , 20141) . 
With respect to quasars, GRBs have many advantages: 
(i) they are already detected at much larger redshift than 
quasars; (ii) they reside in average cosmic regions, less 
affected by local ionization effects or clustering, and 
(iii) their spectrum follows a power-law making con¬ 
tinuum determination much easier. The main limita¬ 
tion of this method is repr esented by the presenc e of 
the host galaxy absorption (IMcOuinn et all l2008h . al¬ 
though t his can be in princ iple separated being less ex¬ 
tended. Chen et all ( 20071) found that ~ 20 - 30% of 
GRBs have small enough hydrogen column density to 
allow a direct measure of absorption from a partially 
neutral IGM. Moreover, the typical Nhi value for GRB 
hosts is expected to de crease with increasing redshift 
(INagamine et alil2008h . This idea is now supported by 
the fact that three out of the four z > 6 GRBs for which 
the intrinsic HI column density has been measured, 
showed Nhi < 10^° cm“^ (see Table 1). In particular, 
a very low Nhi i s foun d in the case of GRB 140515A 
( Chornock et al. , 20141) suggesting that a statistical sam¬ 
ple of bursts with low intrinsic column densities can be 
gathered in future allowing to firmly constrain the reion¬ 
ization process. 

Another effective method to constrain the reioniza¬ 
tion history is based on the statistics of peaks/gaps 
in the afterglow spectrum between Lycr and Ly^S, i.e. 
correspon ding to transmission/absorp tio n regions along 
the LOS ( Paschos & Norm^ 200-5h . 


Gallerani et al 
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(120081) have shown that the distribution of the dimen¬ 
sion of the largest dark gap is sensitive to the assumed 
reionization history provided that a statistical sample of 
~ 20 LOSs is gathered. 

Finally, the study of the HI forest over-imposed on the 
GRB radio afterglow with SKA of a very high redshift 
burst could offer, in princ iple, a powerful tool to study 


the re ionization process (iXu et al.L 1201 It ICiardi et al. 
2015h . However, it requires the detection of an ex 


tremely bright radio afterglow at the level of ~ 10 mJy 
at z > 7 - 8, that may be possible in the case of PopIII 
GRBs. 

Beside constraining the reionization history and its 
variance along different LOSs, GRBs can provide an 
estimate of the typical escape fraction of UV photons 
from early galaxies. The measure of this crucial quan¬ 
tity from the observation of the ecaping Lyman con- 
tinuum radiation is a lready very difficult at z ~ 3 
( Vanzella et all 2012), and may be impossible for the 
small, faint galaxies responsible for the bulk of star for¬ 
mation at z > 6. A statistical sample of GRBs provide 
an alternative way to infer it from the distribution of Ahi 
over many LOSs. Useful constr ains have so far obta ined 
at z = 2 - 4 with /esc < 7.5% ( Fvnbo et al. . 2009h . be¬ 
ing higher redshift studies prevented by the small statis¬ 
tic of the sample. The measure of the escape fraction 
from typically star-forming galaxies during the epoch of 
reionization will provide the missing piece of the puzzle 
in our understanding of how galaxies reionize the IGM. 

4.3. Studying high-z galaxies 

Theoretical models and observational evidence sug¬ 
gests that most of the UV photons that reionize the Uni¬ 
verse are emitted by low mass, faint galaxies missed 
even in the deepest observations by HST. The detec¬ 
tion of these objects is one of the main scientific goal of 
the next generation of space telescopes. However, lit¬ 
tle or no information about the physical properties (e.g. 
metal and dust content) of these objects will be acces¬ 
sible even with JWST. High-z GRBs can provide a use¬ 
ful, complementary tool to investigate the building up 
of metals at these early stages of galaxy formation and, 
in particular, in those objects that provide th e bulk of 
the ionizing photons ( Salvaterra et all 2013 ). Indeed, 
GRB events are produce by the same massive stars that 
emit the FUV photons able to ionize hydrogen. Metal 
absorption line observations of a large sample of good 
signal-to-noise GRB spectra will allow to recover the 
cosmic metal enrichment history, to extend the mass 


metal l icity and FMR to hig her redshift (ISalvaterra et al 


to search for peculiar nucleo-synthesis pattern (see next 
Section). 

Metal absorption lines over-imposed to the afterglow 
radiation have been already detected even z > 6 showing 
that high-z galaxies are already enriched to a few per¬ 
cent solar. In particular, the spect rum of GRB 130606A 
obtained with VLT/X-shooter ( Hartoog et al.L 2014 ) 
shows a great variety of metal elements allowing ac¬ 
curate abundance and gas kinematic studies of the in¬ 
ner region of a z ~ 6 galaxy. Even better and richer 
datasets will be gathered when the 30-m class tele¬ 
scope like E-ELT will become operative. These stud¬ 
ies can be complemented by X-ray spectroscopic ob¬ 
servations of bright GRBs with the next g eneration of 
X-ray instruments ( Camnana et al.[ 2011). I n part ic¬ 


ular, with the Athena satellite (INandra et al.L 12013h it 


will become possible to directly measure the abundance 
patterns in X-ray afterglows allowing, in principle, to 
discriminate betwee n different nucleo-synthesis sources 
( .Tonker et al. . 2013 ). 

Beside metal absorption lines, the observation of H 2 
molecular absorption and of the local dust law can pro¬ 
vide further details about the host enrichment. As dis¬ 
cussed in Section l272l evidence of dust extinction in four 
z > 6 GRBs has been derived by studying the X-ray-to- 
optical SED. Eurthermore, in the case of GRB 140515A 
dust deplitio n is suggested by the peculiar metal abun¬ 
dance ratio ( Hartoog et al.L 2014 ). These observations 
show that GRBs can provide fundamental clues about 
the presence of dust in high-z galaxies. Moreover, they 
can be use to constrain the dust formation channels at 
work at those early epoch. I t is interesting to note 
that for GRB 071025 at z ~ 5 dPerlev et al.Ll2010l) and 


tentat ively for GRB 050904 at z = 6.3 (IStratta et al 


201 Ih evidences for SN synthesized dust has been re¬ 
ported. Tracking the dust enrichment history of first 
galaxies is also very important for better understand the 
PopIII/PopII transition. Indeed, it has been shown that 
even a little amount of dust in the medium will induce 


the fo rmation of low mass stars (iSchneider et al.L 12003 

2 OI 2 I) . 


2013 : Laskaret al.l 2011) , to perform ste l lar po pula- 
tion studies ( Grieco et al. . 2014 : Ma et al. . 20L5h . and 


4.4. Indirect search for PopIII stars 

PopIII GRBs (if any) are expected to be extremely 
rare and difficult to distinguish from other GRB pop¬ 
ulations (see Section I3.2l i. An indirect, but fruitful, 
way to search for the elusive PopIII stars is repre¬ 
sented by the study of the metal and dust compo¬ 
sition of the ISM of distant galaxies enlightened by 
PopI I GRBs. Indeed, peculi ar metal abundance ra¬ 


tios (iHeger & Wooslevi l2002h and dust composition 
(ISchneider et al.L 120041) are foreseen as the result of the 
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enrichment of first massive SN explosions. The detec¬ 
tion of such a signature in the optical-NIR afterglow of 
a very high-z GRBs will provide a strong evidence for 
th e existence of very massive (> 100 Mq) PopIII stars. 


Wang et al.l (120121) found that PopIII enriched gas 


will result in much larger equivalent widths of metal ab- 
sorption lines th an typical PopII enrichment. Moreover, 
Ma et al. ( 2015 ) showed that in principle the [C/O] and 
[Si/O] alone could be enough to distinguish a PopIII en¬ 
riched envirorment, although, in practice, the detection 
of more elements (such as S and Fe) is needed. PopII 
GRBs exploding in a PopIII enrichment medium are ex¬ 
pected to be rare atz = 6(< l%of the GRB popu¬ 
lation at that redshift) but the probability to find them 
increa ses with z being ~ 10% at z = 10 (IMa et al 


2015h . They should reside in small host galaxies with 
M* ~ 10^“® Mq and Z < 10“^ ^ Zq. Thus, the detection 
of a z > 10 burst with a metallicity below this threshold 
will be a strong candidate to be either a PopIII event or 
a PopII GRBs blowing up in a PopIII enriched medium. 


4.5. Others 


For the sake of comp leteness (see iMcOuinn et al. 
2009; Amati et al. . 201 .31 for a more complete list) I 
would like to briefly mention here some other possible 
uses of GRBs in the study of early Uni verse: _ 


c onstrain the dark matte r particle (iMesinger et al. 


2005; de Souza et all 2013h . The redshift distribution 
of GRBs can set limits on the dark matter particle mass, 
mx. Using a sub-sample of z > 4 GRBs mx > 1.6 - 1.8 
keV at 95% confidenc e level is found. 

- non-Gaussianity ( Maio et HI 2012 ). Deviations 

from Gaussianity of the primordial density field will 
translate in different rates of GRBs z » 6. A single 
GRB detection at z > 15 will favor non-Gaussian sce¬ 
narios with a positiv e non-linear parameter. _ 

- radiation field ( Inoue et al. . 2010l : Kakuwa et al. 


20121). The ultraviolet intergalactic radiation field be¬ 
low the Lyman edge energy can cause attenuation in the 
y-ray spectra of GRBs. This may be observable at high- 
z by the Cherenkov Telescope Array (CTA ). 

- magnetic fields (iTakahashi et al.LUOl ih . Pair echoes 
from luminous high-z GRBs may be detectable with 
CTA providing a unique probe of weak intergalactic 
magnetic fields at early stages of structure formation. 


5. Present and future of the fields 


The detection of high-z GRBs is one of the main 
goals driving the design of the next GRB missions 
(I Amati et all l2013h . As already discussed, GRBs are 



10 -« 10 -“ 0.0001 0.001 0.01 0.1 1 


PjE [Ph S-' cm-®] 

Figure 4: Required sensitivity, in terms of minimum peak flux P^e 
that can be detected in the given energy band AE, and field-of-view to 
detect 10 GRB yr“* at z > 8. Different lines conespond to different 
ene rgy bands as labeled in the plot, i.e. to different mission concept. 
See lGhirlanda et all 1^1^ for the details of the calculation. 


very effective in selecting high-z objects and sufficiently 
bright to be detectable up to beginning of the star for- 
mati on activity. While the forthcoming SVOM satel- 
lite ( Godetetali 2012 ) or the proposed LOFT mission 
( Amati et aU 2015 ) will surely increase the z > 6 GRB 
sample, it is also clear that in order to full exploit the 
potentiality of GRBs as a probe of the early Univese, 
a much larger sample of well studied, high-z GRBs 
should be collected. 


Due to redshift scaling, in the observer frame the 
prompt emission of high-z GRBs is expected to peak 
at softer energies than lower redshift bursts. However, 
as the sensitivity of current GRB satellites allow us to 
detect only the bright-end of the GRB LF at high-z, the 
mean observed peak energy is < Lp^obs >~ 150 keV 
similar to low-z bursts. Therefore, while a facility oper¬ 
ating in the soft X-rays will be, in general, more efficient 
in selecting high-z GRBs, a much better sensitivity than 
current instruments should also be foreseen. 


The scientific goals outlined in the previous section 
call for a large statistical sample of well studied high-z 
GRBs. To be competitive with other probes, at least a 
hundred of GRBs at z > 6 with several tens of them ly- 
ing at z > 8 should be detected. Ghirlanda et alJ ( 2015 ) 
computed the expected detection rate of high-z GRBs 
by a generic detector with defined energy band and sen¬ 
sitivity by means of a observational tested population 
synthesis model of long GRBs. Following the results of 
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this work, Figure 4 shows the sensitivity and field-of- 
view of a mission able to detect 10 GRBs yr ' at z > 8 
and operating in different energy bands. For instance, by 
adopting the Swift sensitivity of ~ 0.4 ph s"* cm“^ and 
energy band, we can see that the goal is never reached. 
On the other hand, limiting ourselves to the Swift FOV 
of 1.4 sr, a hundred times better sensitivity than Swift is 
needed. A X-ray instrument will instead require a sen¬ 
sitivity of 0.1 ph s ' cm“^ in the 2-5 keV band for the 
same FOV. It is worth to note that the same mission will 
also detect thousands of bursts at lower redshifts (plus 
other transients). In order to effectively select high- 
z candidates, the presence of a 0.5-1 m infrared tele¬ 
scope on-board should be foreseen (together with fast 
repointing capabilities). This instrument will allow not 
only to promptly measure the redshift but also to per¬ 
form low-resolution (R~ 1000) spectroscopic studies 
(e.g. to identify metal absorption lines) when the after¬ 
glow is at the maximum of its brightness. A rapid dis¬ 
semination of best high-z candidates is also fundamental 
in order to trigger follow-up observations by future fa¬ 
cilities operating at different wavelengths: optical-NIR 
(E-ELT, JWST if still operating), radio (SKA), X-ray 
(ATHENA), TeV (CTA). 

In conclusion, the study of the high-z Universe with 
GRBs requires a X-ray detector with unprecedented 
combination of sensitivy and EOV coupled with an in¬ 
frared telescope to select reliable high-z candidates. The 
THESEUS mission recently proposed for the M4 ESA 
call matches all these characteristics. 

6. Conclusions 

In this paper, I have briefly reviewed the status and 
the prospect for the exploration of the high redshift Uni¬ 
verse with GRBs. A few solid considerations can be 
drawn: 

(i) GRBs do exist at very high redshift and can be 
detected and studied with present-day facilities; 

(ii) high-z GRBs are very similar to low- and 
intermediate-z ones; 

(iii) GRBs are an efficient way to select high-z ob¬ 
jects; 

(iv) high-z GRBs are hosted in faint galaxies that are 
missed even in the deepest surveys; 

(v) high-z GRBs have proved to be an independent 
and powerful tool to study the early Universe. 

Although the current sample is limited to a few 
bursts, future dedicated missions can provide a suffi¬ 
ciently large number of GRBs at z > 6 to study the first 
bound structures in a complementary way with respect 
to galaxy and quasar surveys. 
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